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^ ' Abstract. Two bipolar proto-planetary nebulae, IRAS 16594-4656 and IRAS 17 150-3224, have been detected in various molec- 

ular lines, namely CO, ''CO, HCN and CN, and remain undetected in several other species. CO(/=2-l) and CO(y=3-2) 

■ line profiles are compared to new spectra of similar PPN candidates, previously undetected in CO(y=2-l): CPD-53°5736, 
IRAS 1 7 1 06-3046, IRAS 1 7245-395 1 and IRAS 1 744 1 -24 1 1 . CO(7=2- 1 ) maps of IRAS 1 6594-4656 and IRAS 1 7 1 50-3224 show 
that both PPNe have little separation between blue, centre and red-shifted emission, and also that the CO(y=2-l) emission is 
of a similar size to the telescope beam. Fractional abundances of all detected molecules (except CO) are calculated using the 

^ ' results of CO line modelling and a simple photodissociation model. For those species not detected, upper limits are derived. 

, Comparisons between these fractional abundances and those of other PPNe show that IRAS 16594-4656 and IRAS17150-3224 

are quite under-abundant when compared to molecule-rich sources like CRL618, CRL2688 and OH23 1.8+4.2. As a reason 
for this deficit, the difference in circumstellar envelope/torus density between the molecule-rich sources and the molecule-poor 
sources is proposed, and supported by a chemical model which follows the transition of a circumstellar envelope through the 
' AGB phase and into the PPN phase of evolution. The model includes the effects of UV radiation, cosmic rays and also X-rays. 

, Finally, the post- AGB ages of these two objects (200 - 400 yr) are estimated using CN/HCN and HCN/CO ratios and both ages 

■"■^ — , ■ are found to be in agreement with previous figures cited in the literature, IRAS 17 150-3224 being the younger of the two PPNe. 

^p^, Key words. Astrochemistry - Stars: AGB and post- AGB - Stars: carbon -circumstellar matter -Stars: individual: IRAS 16594- 

O ■ 4656 - Stars: individual: IRAS 17150-3224 

■ 1. Introduction such molecular "richness". The reason for this difference is far 
J> from clear. 

The proto-planetary stage of evolution is one of the shortest in 
rS a star's lifetime and this implies a scarcity of observable exam- Many optical and infrared HST images have helped in the 
C3 pies. Following the method of Olivier et al. (2001), the number study of PPNe (e.g., see Su et al. 2001; Hrivnak et al. 2001, 
surface density of proto-planetary nebulae can be estimated to 1999; Kwok et al. 1998, for details on the two PPNe de- 
be (0.41+1,^1) kpc-2. In comparsion, the density of AGB stars is scribed here), although the depth to which many of the objects 
15kpc"2 and of Main Sequence stars is -2x10*' kpc ^ (Olivier discovered have been studied is minimal. Both molecule-rich 
et al. 2001). Hence knowledge about this phase is limited. PPNe and molecule-poor PPNe appear similar in images - the 
Most of the current understanding of PPNe is derived from a three PPNe mentioned above and both IRAS 16594-4656 and 
handful of objects; mainly the carbon-rich sources CRL618 IRAS 17 150-3224 have some degree of bipolarity, and a narrow 
and CRL2688 and the oxygen-rich object OH231.8-H4.2. All waist. The beginnings of this bipolarity are found in the late 
three show strong molecular lines, have axisymmetric struc- AGB phase, where the first signs of asphericity are seen (e.g., 
tures and molecular tori or disks. However, several hundred Kastner & Weintraub 1994). The shaping of the nebula contin- 
PPNe candidates have been identified (see Kwok 1993, and ues under the influence of a fast superwind, according to the 
references therein), some 34 are reasonably well identified as Generalized Interacting Stellar Wind (GISW) model (Balick 
PPNe (Bujaii-abal et al. 2001), and yet only a handful show 1987). Other efi^ects then play apart in the developing morphol- 
ogy (see the review by Balick & Frank 2002), and one which 
may be important is the inertial confinement of the outflowing 
Send offprint requests to: TIM, wind by a circumstellar torus or disk (Calvet & Peimbert 1983). 
e-mail: Tom.Millar@umist.ac.uk The degree of coUimation produced by this disk would depend 
* Figures A.1-A.3 are only available in electronic form via on the mass of the disk, as well as the momentum involved 
http : //www . edpsciences . org in the high-speed outflow, and other, geometrical effects. This 
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Table 1. Positions and LSR velocities for the PPNe under study. 



TR AS No 


wtllCl IICUllC 


J2000 Co-ords. 


vlsr 






[h:m:s] 


[=:':"] 


[kms-1] 


14488-5405 


CPD-53°5736 


14:52:28.7 


-54:17:43 


-10 


16594-4656 




17:03:09.7 


-47:00:28 


-25 


17106-3046 




17:13:51.7 


-30:49:40 





17150-3224 


AFGL6815S 


17:18:20.0 


-32:27:20 


-1-15 


17245-3951 


OH348.8-2.8 


17:28:04.8 


-39:53:44 





17441-2411 


AFGL5385 


17:47:10.3 


-24:12:54 


-1-110 



atures above the atmosphere ranged between 150 and 800 K, 
depending on frequency and elevation. 

Most observations were carried out with the dual beam 
switching method which gives very flat baselines by placing 
the source alternately in two beams. Beam separation was ll'.S 
in azimuth. CO maps of IRAS17150-3224 were taken using 
the position switching method, to minimise interference from 
interstellar CO emission Unes. The OFF position was chosen to 
be (-10', -11') from the ON position. Map spacing in all cases 
was 1 1", which is approximately half a beamwidth at 230 GHz. 
Calibration was performed with the standard chopper-wheel 
technique, and has an uncertainty of approximately 20% (see 
Schoier & Olofsson 2001, for details). J2000 positions of both 
sources are given in Table 1. 

Intensity scales of spectra in this paper are given in 
main-beam brightness temperature, which is the corrected an- 
tenna temperature divided by the main-beam efficiency (?7mb)- 
Beam efficiencies and FWHM beam widths are summarised in 
Table 2. Individual spectra are shown in Appendix A, and hne 
frequencies and integrated intensities are shown in Table 3. 

Integrated intensities of lines which were not clearly de- 
tected are calculated using the following expression, 

/v < 3cr ^ V2 ^2^ J Avres = 3(r(2n)i/2Av,es. (1) 

where cr is the rms noise in the spectrum, Vexp ffie expansion 
velocity, Avres the velocity resolution of the spectrum and n 
the number of channels covering the line width. This gives a 
somewhat optimistic view of the integrated intensity of ffie line, 
and hence any fractional abundance calculation using this value 
will be strictly an upper limit. Values of calculated in this 
manner are given in Table 3 and indicated by a "less than" sign. 

3. Sources 

The two PPNe presented in this paper are objects for which 
there are currently few millimetre-wavelength spectra and little 
or no molecular information. Both show bipolar morphology 
in HST images and also bright CO emission, wiffi interesting 
features. These two sources were selected on the basis of previ- 
ously pubhshed CO spectra (Loup et al. 1990; Hu et al. 1993). 

3.1. IRAS1 6594-4656 

IRAS 16594-4656 is a proto-planetary nebula which lies almost 
in the plane of the Galaxy (b = -3°3, van de Steene & van 
Hoof 2003) and has a mixed chemistry. In its infrared spec- 
trum there are strong and rarely observed features at 12.6 and 
13.4yum, thought to be due to poly cyclic aromatic hydrocar- 
bons (PAHs) with a high degree of hydrogenation, and further 
PAH features in ffie range 3-13yum (Garcfa-Lario et al. 1999). 
An emission feature between 19 and 23 /im, associated wiffi 
C-rich PPNe by Kwok & Hrivnak (1989), also suggests ffiat 
this source is carbon rich. However, ffiere are also indications 
of crystalline silicates (pyroxenes), with weak features around 
34;um (Garcfa-Lario et al. 1999), and a 10 nm silicate feature 
(Olnon et al. 1986). This suggests ffiat until recently this object 



Table 2. SEST FWHM beam widths and efficiencies at selected fre- 
quencies. 



Frequency 


FWHM 


'7mb 


[GHz] 


["] 




86 


57 


0.75 


100 


51 


0.73 


115 


45 


0.70 


147 


34 


0.66 


230 


23 


0.50 


265 


21 


0.42 


345 


15 


0.25 



change in morphology, from something approximately spher- 
ical to something bipolar, or elliptical, occurs very rapidly at 
ffie end of the AGB phase of evolution (e.g., Kwok et al. 1996; 
Schmidt et al. 2002). The actual period of transition is hard to 
quantify, and there are inherent difficulties in estimates by dy- 
namical means (Zijlstra et al. 2001). 

In this paper, two objects which have been imaged by 
ffie HST and show strong CO emission are studied. Both ob- 
jects, IRAS 16594-4656 and IRAS 17 150-3224 are only de- 
tected in a handful of species, including HCN and CO. These 
are compared to other PPNe newly observed in HCN and 
CO - CPD-53°5736, IRAS 17 106-3046, IRAS 17245-3951 and 
IRAS 1744 1-24 11 (Sect. 4). Furthermore, the molecular prop- 
erties of ffiese two PPNe are studied by calculating fractional 
abundances from SEST data, and comparing with similar ob- 
jects (Sects. 6.2.1 & 7). The comparative under-abundance of 
molecules of IRAS 165 94-465 6 and 17150-3224 is discussed 
in Sect. 7, and arising hypotheses are conffimed by means of 
a chemical model (Sect. 6.3). Finally, individual spectra are 
presented in Appendix A, which is only available in electronic 
form. 

2. Observations 

The observations were carried out between 1998 and 2002 with 
the Swedish-ESO Submillimetre Telescope (SEST), situated 
on La Silla, Chile. SIS receivers were used at 0.8, 1.3, 2 and 
3 mm. During this period the SEST operated three acousto- 
optical spectrometers: one high resolution spectrometer, wiffi a 
bandwidth of 86 MHz and a channel separation of 43 kHz, and 
two wideband spectrometers wiffi bandwidths of about 1 GHz 
and a channel separation of 0.7 MHz. Typical system temper- 
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Table 3. Observed lines. 



IRAS 16594-4656 IRAS 17 150-3224 



Molecule 


Transition 


Frequency 




J Tmbdv 




J Tmbdv 






GHz 


K 


Kkms-' 


K 


Kkms"' 


C3S 


7=15-14 


86.708 


<0.01 


<0.03 






SiO 


7=2-1 


86.847 


<0.01 


<0.03 






HN"C 


7=1-0 


87.091 


<0.01 


<0.03 






C4H 


2n3/2, 19/2-17/2 


87.372 


<0.01 


<0.03 






C2H 


Af=l-0 


87.329 


<0.01 


<0.03 






HCN 


7=1-0 


88.632 


0.02 


0.47 


(1.(12 


0.54 


C3N 


Ar=9-8, 7=19/2-17/2 


89.046 


<0.01 


<0.05 






C3N 


^^=9-8, 7=17/2-15/2 


89.064 


<0.01 


<0.05 






HCO+ 


7=1-0 


89.189 


<0.01 


<0.05 






HC'^CCN 


7=10-9 


90.593 


<0.01 


<0.04 






HCC"CN 


7=10-9 


90.602 


<0.01 


<0.04 






HNC 


7=1-0 


90.664 


<0.01 


<0.04 


_ 


_ 


SiS 


7=5-4 


90.772 


<0.01 


<0.04 


_ 


_ 


HC3N 


7=10-9 


90.979 


<0.01 


<0.04 






CH3CN 


5(1)-4(1) 


91.985 


<0.01 


<0.05 






HC"CCN 


7=12-11 


108.711 


<0.01 


<0.06 






HCC"CN 


7=12-11 


108.721 


<0.01 


<0.06 






"CN 


Af=l-0 


108.780 


<0.01 


<0.06 




_ 


C3N 


A'=ll-10, 7=23/2-21/2 


108.834 


<0.01 


<0.06 






C3N 


Af=ll-10, 7=21/2-19/2 


108.853 


<0.01 


<0.06 






SiS 


7=6-5 


108.924 


<0.01 


<0.06 








HC3N 


7=12-11 


109.174 


<0.01 


<0.06 








"CO 


7=1-0 


110.201 


0.04 


0.06 


0.04 


0.75 


CH3CN 


6(1)-5(1) 


110.381 


<0.01 


<0.11 


<0.01 


<0.09 


C4H 


2n3/2 23/2-21/2 


113.266 


<0.01 


<0.09 


<0.01 


<0.06 


C2S 


8(7)-9(8) 


113.410 


<0.01 


<0.09 


<0.01 


<0.06 


CN 


N=l-0 


113.491 


0.05 


3.63 


<0.01 


<0.06 


C4H 


2n3/2 21/2-19/2 


115.217 


<0.07 


<0.76 


<0.04 


<0.45 


CO 


7=1-0 


115.271 


0.61 


1.78 


0.25 


5.05 


SiCz 


5o,5 "4o,4 


115.382 


<0.07 


<0.76 


<0.04 


<0.45 


SiO 


7=3-2' 


130.269 


<0.01 


<0.06 


_ 


_ 


H2CO 


2 12" 1 1.1 


140.840 


<0.01 


<0.07 






SiCj 


62,5 "52,4 


140.920 


<0.01 


<0.07 


_ 


_ 


H"CCCN 


7=16-15 


141.062 


<0.01 


<0.07 






CS 


7=3-2 


146.969 


<0.01 


<0.04 






CH3CN 


8(0)-7(0) 


147.175 


<0.01 


<0.04 


_ 




"CO 


7=2-1 


220.399 


0.10 


1.34 


0.13 


2.91 


CH3CN 


12(0)- 11(0) 


220.747 


<0.01 


<0.16 


<0.01 


<0.17 


CN 


Af=2-1 


226.875 


0.11 


9.91 


<0.01 


<0.10 


CO 


7=2-1 


230.538 


1.75 


32.60 


0.65 


12.84 


CO+ 


7=2-1 


235.790 


<0.01 


<0.12 






SiS 


7=13-12 


235.961 


<0.01 


<0.12 






CH3CN 


13(0)-12(0) 


239.138 


<0.01 


<0.16 






CH3C2H 


14(0)-13(0) 


239.252 


<0.01 


<0.16 






C3H2 


44,1-33,0 


265.759 


<0.02 


<0.36 






HCN 


7=3-2 


265.886 


0.11 


2.85 






H"CN 


7=4-3 


345.340 


<0.04 


<0.80 


<0.06 


< 1.28 


HC3N 


7=38-37 


345.610 


<0.04 


<0.80 


<0.06 


<1.28 


CO 


7=3-2 


345.796 


2.05 


41.20 


0.78 


18.19 



was oxygen-rich, becoming carbon-rich shortly before the tran- 
sition to the post-AGB phase (Garcia-Lario et al. 1999), which 
occurred some 370 years ago (van der Veen et al. 1989). This 
hypothesis was strengthened by the lack of a detection in OH 
(te Lintel Hekkert et al. 1991; Silva et al. 1993) and no SiO 



maser emission (Nyman et al. 1998). In evolved stars with a 
mixed chemistry such as this, the silicates are found in a torus 
and the PAHs in the polar flow (Matsuura et al. 2004). 

Several authors describe HST images of this PPN (Hrivnak 
et al. 1999, 2000; Garci'a-Lario et al. 1999), which show the 
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presence of a bright , B7-type (van de Steene et al. 2000) cen- 
tral star, surrounded by a multiple-axis bipolar nebulosity with 
a complex morphology. The nebula is seen at an intermediate 
inclination, rather than edge-on. The optical size of the nebula 
is given as 5"xH"(Garcia-Lario et al. 1999). 

The only CO spectrum of this source previously published 
is a CO (7=1-0) spectrum from Loup et al. (1990), who de- 
tected broad CO emission corresponding to that from a circum- 
steUar shell expanding with a velocity of 16 km s"^ 

3.2. IRAS1 71 50-3224 

IRAS 17 150-3224 (AFGL68 15S) is a young PPN which left the 
AGB some 150 (Hu et al. 1993) to 210 years ago, and possi- 
bly came from a high-mass progenitor (Meixner et al. 2002). 
Note however, that van der Veen et al. (1989) give a dynamic 
timescale of 800 yr for this object. Structurally, it is a bipolar 
nebula, with a nearly edge-on (82°) dusty torus (Kwok et al. 
1996), expanding at 1 1 km s"' (Weintraub et al. 1998) around a 
G2-type central star. The optical size of the nebula is approxi- 
mately 12" X 9" (Hu et al. 1993; Kwok et al. 1996; Garcia- 
Hernandez et al. 2002). A newly-discovered equatorial loop 
seen in near-infrared, but not optical, images, is discussed by 
Su et al. (2003). A faint halo (AGB mass-loss remnant) can be 
seen in images from Kwok et al. (1998). 

IRAS17150-3224 shows molecular hydrogen emission, but 
not hydrogen recombination, indicating that the degree of ion- 
isation is small (Garcia-Hernandez et al. 2002). The emission 
is Ukely due to shock-excitation, as is often seen in strongly 
bipolar nebulae. This object is similar to CRL2688 (Sahai et al. 
1998) in that it has Hi emission, and a spectral type later than 
A (CRL2688 has a spectral type of F2, Cohen & Kuhi 1977). 
IRAS 17 150-3224 shows OH maser emission (Hu et al. 1994), 
but not SiO maser emission (Nyman et al. 1998). It also has 
been detected in the 3.1 jum water Une (van der Veen et al. 
1989), but not in H2O maser emission (Zuckerman & Lo 1987). 

Previous CO observations are hmited to the paper of Hu 
et al. (1993), who present both CO (7=1-0) and (7=2-1) spec- 
tra. 

4. CO observations 

4.1. IRAS 16594-4656 

The CO (7=2-1) and (7=3-2) line profiles presented in Fig. la 
are affected by interstellar emission in the OFF position (be- 
tween -12 and -5 kms ' approximately). The CO (7=3-2) pro- 
file is taken in dual beamswitch mode whereas the CO (7=2-1) 
profile was taken in position switch mode with an OFF posi- 
tion -30" away in R.A. in order to minimise the eflFect of the in- 
terstellar Une. Despite this, a high signal-to-noise ratio aff'ords 
some indication of structure in this source, shown by features 
in the line profile. A roughly paraboUc profile and a source 
size comparable to the beam shows that the CO emission is 
probably optically thick. An outflowing wind of 14kms~^ is 
indicated by the fine profile, and there seems to be an appre- 
ciable line wing on the blue side, possibly indicating a sec- 
ond wind of up to 25kms"^ The CO (7=2-1) map of this 



source (Fig. 3) shows a reasonably symmetric pattern of emis- 
sion. The red/blue- shifted emission contour plot (Fig.2) does 
not show the red-shifted emission, due to the interstellar inter- 
ference. The blue-shifted emission, taken in the interval [-50,- 
35] km s"\ lies directly on top of the centre emission, taken in 
the interval [-35,-12]kms"', indicating that there is little or no 
separation in emission regions. There is some emission at +20" 
in R.A., and an elliptical Gaussian fit to the integrated intensity 
map shows that the FWHM ellipse has a major axis of ~25" 
and a minor axis of ~22.5" (position angle of 57°), very simi- 
lar to the size of the SEST beam (23") at 230 GHz. Hence the 
source is on the Umit of being resolved. 

4.2. IRAS1 71 50-3224 

The CO (7=2-1) spectrum (Fig. lb) of IRAS 17 150-3224 was 
taken in position switch mode, with an OFF position of -1-45" 
in R.A. The CO (7=3-2) was taken in beamswitch mode. Both 
CO line profiles of IRAS 17 150-3224 shown in Fig. lb are dis- 
tinctly triangular in shape, similar to those shown by 89 Her, 
the Red Rectangle, and M2-9 (and also IRAS 1744 1-241 1, pre- 
sented in this paper). Here the triangular shape is afi'ected by 
an interstellar Une in the blue-shifted wing (between and 
10km s ' approximately). The CO Unes indicate a wind of 
some 30 kms"' coming from the source, and yet a paraboUc 
fit to the Une indicates a wind of 14.5 km s"' . This latter wind 
is ostensibly the AGB wind of the star, with the higher velocity 
emission probably coming from an inner swept-up shell, mov- 
ing at a velocity intermediate to the AGB wind and a post- AGB 
or super-wind. It must be noted, however, that a parabolic fit to 
a triangular lineshape will usually underestimate the expansion 
velocity of the outflowing wind. 

The CO (7=2-1) emission map (Fig. 4) is taken in posi- 
tion switch mode to minimise the effects of the interstellar 
interference, and although hampered slightly by poor point- 
ing, seems to show a reasonably symmetric distribution which 
does not appear to be resolved by the telescope beam. The 
contour map (Fig. 2) does not show any separation in the 
source. The red-shifted emission was taken from the interval 
[30,40] km s~\ the centre emission from [5,30] kms"' and the 
blue-shifted emission from [- 1 8,-6] km s~' . A two-dimensional 
elliptical Gaussian fit to the centre emission gives an extent of 
~20"x26" (position angle of -160°) for this source. Again, this 
is similar to the size of the SEST beam at 230 GHz. 

4.3. Other sources observed in CO 

Spectra of select sources which have been hitherto unobserved 
and/or unpublished in the CO (7=2-1) line are presented in 
Fig. 5. HCN (7=1-0) spectra have also been taken (not pre- 
sented), and HCN/CO intensity ratios are derived (Table 4). 

CPD-53°5736. Although not clear from the noisy spec- 
trum, CPD-53°5736 seems to show the presence of a wind ex- 
panding at 15 km s~' . 

IRASl 7106-3046. This object shows a strange profile: some 
low-level emission, perhaps a wind of ~ 15 kms"', and a thin, 
narrow spike. This narrow spike (of width ~2kms"') could 
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Tmb (K) 




T„b (K) 
2 

1.5 



0.5 



I ' ' ' I 

IRAS1 6594 
CO J=2-1 H 



IRAS1 6594 
CO J=3-2 




Tmb (K) 



-40 -20 
Velocity (km/s) 



0.2 


Tmb (K) 



0.5 



I ' ' ' 

IRAS17150 
CO J=2-1 




IRAS17150 - 
CO J=3-2 




-20 



20 40 

Velocity (km/s) 



60 



Fig. l.a) Top. A position-switched high- 
resolution CO J=2-l spectrum toward 
IRAS 16594-4656. Bottom. A beam- 
switched CO 7=3-2 spectrum toward 
IRAS 16594-4656. 

b) Top. A position- switched high-resolution 
CO 7=2-1 spectrum toward IRAS17150- 
3224. Bottom. A beam-switched CO 7=3-2 
spectrum toward IRAS17150-3224. 



correspond in some way to the circumstellar disk observed by 
Kwok et al. (2000). The narrow feature seems to be offset from 
the centre of the broad feature, although the signal-to-noise ra- 
tio is low. 

IRAS17245-3951 .Low signal-to-noise means that this pro- 
file is hard to quaUfy. The profile suggests an expansion veloc- 
ity of 15 km s"' , but this value does have a considerable uncer- 
tainty. 

IRAS17441-241 1 . This source has a triangular lineshape, 
similar to IRAS 17 150-3224. It is somewhat asymmetric, and 
quite narrow in velocity (20kms~' in width). There are no ap- 
preciable hue wings. 

4.3.1. Ihcn/Ico ratios 

The Ihcn/Ico ratios shown in Table 4 (~ 0.01 - 0.05) look to 
be reasonably consistent throughout the sample of PPN can- 
didates. To put these figures into context, a sample of oxygen 
stars selected by Lindqvist et al. (1988) has an average Ihcn/Ico 
ratio of 0. 124, whereas the carbon star sample of Olofsson et al. 
(1990) has an average Ihcn/Ico ratio of 0.623. Both these sam- 
ples were made up of AGB stars; the values for PPNe seem to 
be smaller. IRAS07134H-1005 and IRAS 19500- 1709, both C- 
rich PPNe, have Ihcn/Ico ratios of 0.15 and 0.08, respectively 
(Bujarrabal et al. 1992). Evolved PNe generally have IhcnAco 
ratios of less than 0.1 (c.f., Bachiller et al. 1997a). 



5. CO line modelling 

A detailed, non-LTE, radiative transfer code based on the 
Monte Carlo method was used in order to model the observed 



molecular line emission. The code assumes the emitting en- 
velope to be spherically synometric and expanding at a con- 
stant velocity. The thermal balance equation, where CO is as- 
sumed to be the main molecular coolant, is solved together with 
the molecular excitation in order to obtain the kinetic temper- 
ature structure. The CO observations are used as constraints 
in the modelling to determine also the density and velocity 
fields through a ;jf'^-analysis. This approach was also adopted 
by Woods et al. (2003b) when modelling a sample of high 
mass-loss rate carbon stars. More details on the radiative trans- 
fer modelling can be found in Schoier & Olofsson (2001). The 
code has been benchmarked to a high accuracy (van ZadelhofF 
et al. 2002). 

In order to derive the mass loss rate which, together with the 
expansion velocity, sets the density scale in the wind through 
the continuity equation, a CO abundance relative to H2 must be 
adopted. For IRAS 16594-4656 a value of 1 x 10 ^ which is 
typical of a C-rich AGB envelope, is used, whereas for the O- 
rich object IRAS17150-3224 a value of 2 x lO"'' is adopted. 
The derived envelope properties such as the mass loss rate 
(M), expansion velocity (Vexp), extent of the CO envelope (/?e) 
and ^^CO/^^CO-ratio are reported in Table 5 together with the 
adopted distance (D), luminosity (L) and initial CO fractional 
abundance (/co). Also shown is the /z-parameter that deter- 
mines the amount of heating in the envelope due to momen- 
tum transfer due to dust-gas collisions (see Schoier & Olofsson 
2001, for details). 

In Fig. 6 the best fit model for IRAS16594-4656, using a 
mass loss rate of 1 X 10"^ Moyr"' and a '^CO/'^CO-ratio of 
30, is overlayed on the observed spectra. The total integrated 
intensities in the lines are well reproduced as are the overall 
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Fig. 2. Position-position contour emission maps of IRAS 16594-4656, and IRAS 17 150-3224. Blue contours show blue-shifted emission, black 
contours emission from the centre of the line profile, and red contours red-shifted emission. 

Table 4. Observed lines. 
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J T'mbdv 
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Kkms-' 


HCN(y=l-0) 

co(y=2-i) 

Ihcn/Ico 


0.09 


<0.07 
2.07 
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0.02 
1.75 


0.47 
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O.OI 


0.30 


1.84 


0.02 
0.65 


0.54 
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0.04 


0.04 


<0.09 
0.58 
<0.15 


0.05 
0.57 


0.32 
6.24 
0.05 



line profiles. There is confusion due to interstellar lines at the 
red-shifted edge of the line profiles. Also, on the blue-shifted 
edge there are signs of a second weak component, possibly a 
higher velocity wind due to the present-day mass loss. 

A good fit to CO observations is also obtained for 
IRAS 17 150-3224 using a mass loss rate of 3 x 10"^ Moyr"' 
(Fig. 7). The '^CO/'^CO-ratio is 7, significantly lower than for 
IRAS 16594-4656 but consistent with IRAS 17 150-3224 being 
O-rich (Abia et al. 2001). The '^CO spectra for this source 
show signs of excess emission at larger expansion veloci- 
ties possibly indicating the presence of a faster moving wind. 
However, the blue side of the emission is confused by interstel- 
lar line contamination. 

The /z-parameter of 10 derived for IRAS 16594-4656 is sig- 
nificantly larger than that of IRAS 1 7 1 50-3224 and that which is 
typically derived for high mass-loss rate AGB-stars, h ~ I - 2, 
(Schoier & Olofsson 2001; Woods et al. 2003b). The kinetic 
temperature derived through the envelope of IRAS 16594-4656 
is about a factor of two, on average, larger than that obtained for 



the carbon star IRC-i-10216 (Schoier & Olofsson 2000, 2001). 
Such a large discrepancy is hard to explain due to, e.g., uncer- 
tainties in the adopted distance, and could instead indicate that 
the properties of the dust grains are different in IRAS 16594- 
4656 or that there is an additional source of heating of the enve- 
lope. A further possibility is that the spectra are contaminated 
by a second, warmer component such as the bipolar present day 
mass-loss seen in optical images at smaller spatial scales. 

The winds of these PPNe are assumed to resemble those of 
their progenitor AGB-stars, i.e., spherical symmetry and a con- 
stant expansion velocity are assumed. Given that both sources 
left the AGE ~ 200 -400 yr ago (see Sect. 7.3) this assumption 
should be valid for distances > 1 x 10'^ cm, since the material 
outside this radius would once have made up the AGB wind. 
However, there is the possibility that a high velocity bipolar 
wind is penetrating this remnant AGB shell. The CO emission 
from the lower rotational transitions observed here are formed 
mostly at distances larger than 5 x 10'^ cm and thus mainly 
probe the outer, AGB-type, wind. For other molecular species 
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Fig. 3. A map of the CO 7=2-1 line around IRAS 16954-4656. The beam-spacing is 11". The spectra are affected by interstellar line contami- 
nation in the red-shifted line wing. 



(e.g., HCN) contributions from the present-day wind might be- 
come significant. 



However, a more detailed treatment of the wind character- 
istics will require high-resolution, multi-transition observations 
that will become a reahty with the upcoming APEX/ ALMA 
telescopes. Deviations from the AGB wind model will become 
more apparent in the high-/ CO hues, both in the form of the 
line profile (broader wings) and in the intensity. The AGB wind 
model predicts that the higher-/ CO lines (7=4-3 to 13-12) 
should all have an integrated intensity of about 30Kkms ' 
and 12Kkms-i for IRAS 165 94-465 6 and IRAS 17 150-3224, 
respectively, when observed with the APEX telescope. Also, 
the effect of dust emission in the excitation of the molecules, 
which is not taken into account in the present analysis, needs 
to be investigated. 



6. Molecular evolution 

6.1. Molecules 

The molecules seen in PPNe give information on various pro- 
cesses which have occurred in the history of the particular PPN. 
HCN and CN can be used to give an idea of the evolutionary 
status of a PPN: in AGB stars, HCN is usually more abundant 
than CN. As HCN becomes increasingly photodissociated, its 
daughter, CN, becomes increasingly abundant; i.e., one would 
expect the CN/HCN ratio to increase with evolutionary age, in 
the early post-AGB phase (see Fig. 8). A similar effect is seen 
in the comparison of HCN and HNC: In CRL618, for exam- 
ple, the column densities of HCN and HNC are comparable 
(Herpin & Cemicharo 2000), whereas in IRC-i- 10216 (a less 
evolved object) HCN is far more abundant that HNC. In both 
these objects HNC is formed via dissociative recombination of 
HCNH-^. 

The HCN/CO ratio should also give some indication of age, 
since photodissociation would cause this ratio to decrease with 
time. However, it may be difficult to compare a number of ob- 
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Fig. 4. A map of the CO .1=2-1 line around 1RAS17150-3224, taken in beamswitch mode. The beam-spacing is 11". The spectra are affected 
by interstellar line contamination in the blue-shifted line wing. 

Table 5. Adopted and derived modelling parameters of the two PPNe. 



IRAS No. 


D 


L 


fco 


M 






h 


'^CO/"CO 




[kpc] 


[Lo] 




[Moyr-'] 


[kms-'] 


[cm] 






16594-4656 


1.80" 


4900 


1x10"^ 


1x10-5 


14.0 


2.9x10" 


10.0 


30 


17150-3224 


2.42* 


11000 


2x10-^* 


3x10-5 


14.5 


2.4x101'' 


1.8 


7 



(a) van der Veen et al. (1989). 

(b) Bujarrabal et al. (2001) 



jects which includes a mixture of O-rich and C-rich stars, since 
they may have had a different HCN/CO ratio when they started 
to move off the AGB. 

The presence of CS can be used in a similar way. CS is 
abundant in AGB envelopes, and although it can be enhanced 
by shocks in PPNe (Kasuga et al. 1997; Mitchell 1984) is gen- 
erally of lower abundance, and unobserved in PNe until re- 
cently (Woods & Nyman 2004). HC3N, SiO and SiCa demon- 
strate the same sort of behaviour, decreasing in abundance go- 
ing into the PPNe phase. 



HCO^ becomes greatly enhanced in the PPNe phase, and 
yet it is very difficult to explain why. Theoretical models of 
post-AGB objects all fail to produce enough HCO"^ by at least 
an order of magnitude, if not more (e.g., Howe et al. 1994; 
Ali et al. 2001; Woods et al. 2003a). HCO+ is only efficiently 
formed by ionising radiation, and yet the levels of cosmic 
ray ionising radiation required to produce the desired levels 
of HCO^ result in the destruction of most other molecules. 
This is due to the fact that cosmic rays induce UV photons in 
their interaction with H2; these photons dissociate and ionise 
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Fig. 7. Best fit CO model of IRAS 17 150-3224 (solid line) overlayed on observations (histogram) using a mass loss rate of 3 x 10 Mq yr ' and 
a '^CO/^CO-ratio of 7. 



molecules. Woods et al. (2003a) considered whether X-rays 
might be the responsible for the high level of HCO^. Their 
treatment was quite general - a treatment where X-rays are 
solely responsible for ionisation (and not dissociation) might 
prove interesting. This approach was taken in the chemical 
model explained in Sect. 6.3. However, a more thorough treat- 
ment needs to be made. 

HCO^ can also form in photon-dominated regions (PDRs), 
in warm gas (800 K), through reactions involving CO^, which 
itself is formed through reactions between and OH 
(Hasegawa et al. 2000). Moreover, shocked gas can be iden- 
tified by the presence of HCO^ (as well as SiO) (e.g., Mitchell 
& Watt 1985), and although high levels of HCO+ do not nec- 
essarily imply shocks, the presence of high levels of SiO does 
give a much better indication of the presence of shocked gas. 

The '^CO/'^^CO is low (< 10) in several PPNe (e.g., 
Kahane et al. 1992) compared to 10 < '^CO/'^^CO < 25 in 
evolved PNe (Palla et al. 2000; Bachiller et al. 1997a), 
10 < '^CO/'^CO < 40 in compact PNe (Josselin & Bachiller 
2003) and 10 < '^CO/'^CO < 50 in AGE stars (e.g., Kahane 
et al. 1988; Greaves & Holland 1997; Schoier & Olofsson 
2000). From our radiative transfer analysis we derive 
i^CO/i^CO ratios of 30 and 7 for the PPNe IRAS 16594-4656 
and IRAS 17 150-3224, respectively. This ratio may be an indi- 
cator of evolution (since '^CO self-shields less than '^CO), but 
it is more likely to be an indicator of a star's nucleosynthesis 
history which may be very different for otherwise apparently 
similar objects (e.g., Kahane et al. 1992, 1988). 



6.2. Fractional abundances 

6.2.1 . Calculation of fractional abundances 

In Sect. 5 the mass loss rates, and hence the radial distribu- 
tion of the number density of H2, for the sample PPNe were 
determined from a detailed radiative transfer analysis of the 
observed CO line emission adopting a particular CO/H2 abun- 
dance ratio. Here, fractional abundances for the remaining 
molecular species observed were calculated assuming optically 
thin emission. The radiative transfer for these lines is treated in 
the same way as described in Woods et al. (2003b), and the re- 
sults are shown in Table 6. A constant excitation temperature of 
25 K is assumed throughout the emitting region. Upper limits to 
fractional abundances are given in Table 7, and are calculated 
using the integrated intensity from Eq. 1 . The same molecular 
data as used in Woods et al. (2003b) are adopted here. 

Inner and outer radii of molecular distributions are deter- 
mined in much the same way as Woods et al. (2003b), using a 
photodissociation model. Here, however, there is not the bene- 
fit of good interferometric data. Hence simplifying assumptions 
have to be made about certain molecular distributions. The dis- 
tribution of all parent species, with the exception of CO, is 
calculated from the photodissociation model (observed parent 
species are HCN, CS, SiO, SiS). The distribution of CO comes 
from the radiative transfer modelling, as described in Sect. 5. 
Self-shielding of CO is taken into account, in the manner of 
Mamon et al. (1988). All ionisation products or products of the 
circumstellar chemistry are assumed to take on the distribution 
of C2H, with the exception of CN, which has a slightly more ex- 
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Table 7. Calculated upper limits to fractional abundances, with respect 
to m(H2). 



IRAS 16594-^656 IRAS17150-3224 



Molecule /x /x 

CO+ 1.1 (-8) 

HCO* 1.2 (-7) — 

H2CO 1.1 (-6) — 

CS 6.5 (-8) — 

C2S 3.4 (-7) 9.5 (-8) 

C3S 4.1 (-7) — 

SiO 5.1 (-8) — 

SiS 3.2 (-7) — 

SiC2 4.7 (-7) 2.1 (-6) 

CN Det." 1.1 (-8) 

"CN 5.0 (-7) — 

H"CN 2.4 (-8) 1.3 (-8) 

HNC 1-8 (-7) — 

HN"C 1.9 (-7) — 

HC3N 2.7 (-7) — 

H(i3)C3N* 2.7 (-7) — 

H"CCCN 2.0 (-7) — 

C3N 4.2 (-7) — 

CH3CN 8.2 (-8) 5.3 (-8) 

C2H 3.1 (-6) — 

CH3C2H 3.5 (-6) — 

C3H2 2.5 (-7) — 

C4H 4.5 (-6) 1.4 (-6) 



(a) signifies a detection has been made of that particular line in that 
particular source (see Table 6). 

(b) signifies a blend of HCC"CN and HC"CCN. 



pansive distribution. This should not have too great an effect on 
results since there is only a linear dependence on the difference 
between inner and outer radius. These are very straightforward 
assumptions, which ignore effects of circumstellar chemistry, 
and freeze-out onto grains, for example. Again, isotopomers 
are assumed to have the same distribution as their more abun- 
dant forms. 

6.2.2. Uncertainties 

As discussed in Woods et al. (2003b), there are inherent uncer- 
tainties in the approach taken here. Typical distance estimates 
can vary by up to a factor 2, and this influences the mass-loss 
rate derived from the radiative transfer modelling (see Sect. 5). 
When combined with the uncertainties in the choice of inner 
and outer radii for the molecular distributions, the overall error 
in fractional abundance varies with ~D~^~^, where D is the dis- 
tance. Given also that there are likely to be errors introduced 
by the choice of excitation temperature, an overall error of a 
factor 5 or so is to be expected in fractional abundance esti- 
mates. Moreover, some species other than CO are expected to 
be optically thick, most notably HCN, and the optically thin ap- 
proximation would give systematically too low fractional abun- 
dances in those cases. As mentioned in Sect. 5, contribution to 
the line intensities from the present-day, non- spherically sym- 
metric, high-velocity, mass loss might become significant for 



some molecules. A full treatment of this problem, including 
the effects of radiative excitation due to dust emission, will re- 
quire high spatial resolution observations of a large number 
of molecules. With these caveats, molecular fractional abun- 
dances are believed to be order of magnitude estimates. 

6.3. The chemical model 

To investigate the changes in fractional abundance of species 
during the late AGB and PPN phases, a chemical model of 
the circumstellar envelope was used. Results of this model 
are shown in Fig. 8. The chemical model of IRC+10216 con- 
structed by Millar et al. (2000, MHB) has been adapted to in- 
vestigate the chemistry, particularly the photochemistry, as an 
AGB star finishes its phase of mass loss, and moves into the 
PPN phase. Several simplifying assumptions are made, and the 
MHB model has been changed to include X-rays, and a mini- 
mal X-ray chemistry. 

6.3.1. Model parameters and simplifying assumptions 

Several parcels of gas and dust in the circumstellar envelope 
(CSE) of an AGB star are assumed to flow outward at a uni- 
form velocity of 14kms from initial radial distances of 1, 
4 and 7xl0'^ 1, 4 and 7x10'*" and 1 and 4 xlO'^cm from 
the central star. A mass-loss rate of 1x10"^ MOyr"' is chosen, 
in accordance with the molecular line modelling in Table 5. 
A spherical geometry implies that the density of the homoge- 
neous outflow drops as r^^. Mass-loss is assumed to end with 
the start of the model. 

Envelope heating. During the expansion, the "kinetic tem- 
perature" of the parcels varies according to r=ro(ro/r)'^-^^, 
where Tq increases from 100 K to 1 000 K over the first 1 000 yr 
of the model (chosen to represent the PPN phase). This heat- 
ing is due to the increasing temperature of the central star, as it 
evolves towards a white dwarf. 

Stellar UV field. Following the increase in stellar temper- 
ature, the intensity of the stellar UV field also increases. An 
ionisation rate similar to that of the interstellar UV field is as- 
sumed at the start of the model. At the end of the first 1 000 yr 
this value has increased by a factor of 2x10^ (c.f., Herpin & 
Cemicharo 2000). 

X-rays. X-rays from the central star or companion (Woods 
et al. 2002) are treated in a very simple way in this model. The 
X-ray extinction law of Deguchi et al. (1990) is used (see their 
Eq. B3), and an X-ray energy of 1.2keV is assumed (/I ~1 nm). 
X-rays only play an ionising role in the chemical reaction net- 
work, and these reactions are given the same rate coefficients 
as the equivalent cosmic -ray reactions. A very low ionisation 
rate of 10 s ' is chosen initially; this has risen to a rate of 
10"^^ s~^ by the end of the model run. 

6.3.2. Chemical reaction network 

The reaction network is the same as that used in Woods et al. 
(2003a), with the addition of a small number of X-ray ioni- 
sation reactions. At present the network includes 3897 reac- 
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Table 6. Calculated fractional abundances, with respect to «(H2). 
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(a) "CO fractional abundances derived from full radiative transfer (Sect. 5). 

(b) signifies a non-detection, but the calculation of an upper limit to the fractional abundance has been made (see Table 7). 



Table 8. Comparison with fractional abundances in other PPNe, arranged in order of decreasing stellar temperature. 
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(a) This work, (b) Bujarrabal et al. (1994) (c) Bachiller et al. (1997b) (d) Bujarrabal et al. (1988) (e) Sopka et al. (1989) (f) Fuente et al. 
(1998) (g) Nguyen-Q-Rieu & Bieging (1990) (h) Kasuga et al. (1997) (j) Sahai et al. (1994) (k) Using the results of Bujarrabal et al. 
(1992, 2001). (1) Morris et al. (1987) (m) Sanchez Contreras et al. (1997) (n) Omont et al. (1993) (o) Lindqvist et al. (1992) (p) Nercessian 
et al. (1989) (q) Using the results of Olofsson & Nyman (1999). (r) Woods et al. (2003b) 



lions among 407 species in 6 elements. Initial fractional abun- 
dances are determined via the use of an AGB steady-state 
model, with a mass-loss rate of lxlO~^ M0 yr ' and expansion 
velocity of 14kms '. For example, a parcel starting at a ra- 
dius of IxlO'* cm in the model presented here will have initial 
abundances which are the same as a parcel of "AGB" material 
at the same radial distance. The initial abundances of parent 
molecules for the AGB model are chosen to represent a carbon- 
rich environment, and values are taken from MHB. 

7. Discussion 

7. 1. Comparison of observations witli cliemical model 

The chemical model generally shows good agreement with the 
fractional abundances derived here. Given the order of mag- 
nitude error margins (Sect. 6.2.2), both HCN and CN reach 
the calculated fractional abundances for IRAS 16594-4656 dur- 
ing the first 1 000 yr of the model (i.e., the PPN phase). Upper 
limits for other observed species (e.g., HNC, HC3N, C2H and 



C4H) are also in agreement with the model in this region. CS 
is predicted to be more abundant than observed. The model 
also shows that a few species (e.g., CO, C2H and CN) survive 
at observable levels after, say, 1 000 yr of post- AGB evolution 
(Fig. 8). 

Perhaps the most interesting result of the model is the vari- 
ation of the CN/HCN, HNC/HCN and HCN/CO ratios. As ex- 
pected, the CN/HCN ratio is seen to rise during the early post- 
AGB evolution (taking a mean through the envelope. Fig. 8). 
Similarly, the HNC/HCN ratio also increases. The HCN/CO 
ratio declines significantly with post-AGB age. When these 
three behaviours are combined, a reasonably accurate age can 
be given. This approach is taken later, in Sect. 7.3. 

7.2. Comparison of observations witf) otiier 

observations 

IRAS 16594-4656 is very under-abundant in molecules in com- 
parison to C-rich AGB stars. Generally, fractional abundances 
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Fig. 8.Results of the chemical model of 
IRAS 16594-4656. The figure shows the 
variation of the fractional abundances (or ra- 
tios, in the lower row of axes) of selected 
species radially through the expanding cir- 
cumstellar envelope at different points in 
time (measured in years since the cessation 
of mass-loss). The shading marks the pas- 
sage of the inner edge of the circumstel- 
lar envelope. As can be seen at long times 
(5000yr), only CO is reasonably abun- 
dant. The three sets of axes labelled "Oyr" 
show the results of a standard AGB model 
(see 6.3.2). 



derived here, and presented in Table 6 (and the upper lim- 
its shown in Table 7), are an order of magnitude deficient 
when compared to C-rich AGB stars, such as those surveyed 
in Bujarrabal et al. (1994) and Woods et al. (2003b). Two ex- 
ceptions to this are CN and '^CO, which show good agreement 
with the fractional abundances derived for IRC+ 10216. 

The lack of detection of SiO would seem to confirm that 
IRAS16954-4656 is C-rich, and that any silicate grains present 



(see Garcia-Lario et al. 1999) are not being destroyed by 
shocks because either shocks might not be present, or the 
shocks which are present are not strong enough (c.f., van de 
Steene & van Hoof 2003). Another explanation may be that the 
shocked regions are very small, in which case high abundances 
of SiO would be locaUsed, and SiO emission diluted by the 
beam. 
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Conversely, IRAS 17 150-3224 seems (chemically) very 
much like an O-rich AGB star, such as those observed by 
Bujarrabal et al. (1994), or Lindqvist et al. (1988), for exam- 
ple. Although the fractional abundances of only two molecules, 
i^CO and HCN, in IRAS17150-3224 are suitable for compari- 
son, they are very similar to the O-rich AGB stars in both AGB- 
star samples. 

So, IRAS16594-4656 and IRAS17150-3224 seem to have 
different chemical properties from each other, and both seem 
far from the well-known molecule-rich PPNe, like CRL618, 
CRL2688 and OH23 1.8-1-4.2. IRAS16594-4656, supposedly 
a C-rich object, is poor in molecules and, as can be seen 
from Table 8, is at least an order of magnitude less abun- 
dant than CRL618 and CRL2688 in many molecules. However, 
IRAS 16594-4656 has generally higher fractional abundances 
than Ml-16, which is widely accepted to be a standard PN. 
An exception to this is CN - a fractional abundance of CN 
similar to those in CRL618 and CRL2688 may indicate that 
IRAS 16594-4656 is a somewhat evolved PPN, but not quite 
part of the PN regime yet. The theoretical work of Woods et al. 
(2003a) shows that fractional abundances of CN can remain 
high for a period of ~ 100 years in the PPNe phase, whilst other 
molecules are destroyed, if molecular material is to be found in 
the shielded environment of a circumstellar torus. 

Bachiller et al. (1997a) present a homogeneous sample 
of molecular abundances in post-AGB objects; however, both 
IRAS 16594-4656 and IRAS 17 150-3224 are too low in molec- 
ular abundance to present much of a comparison. These two 
PPNe are closest in fractional abundances to the younger ob- 
jects in the Bachiller et al. (1997a) sample. 

Using the fractional abundance of HCN derived from the 
7=1-0 line, since this line is least likely of the two observed 
to be optically thick, a CN/HCN ratio of < 10 is derived for 
IRAS 16594-4656. This value is quite uncertain, and dependent 
very much on the factors discussed in Sect. 6.2.2. A full ra- 
diative transfer treatment of the HCN(./=l-0) line does in fact 
predict a fractional abundance 2-3 times larger than that de- 
rived via the optically thin assumption (Sect. 6.2.2), imply- 
ing that IRAS 16594-4656 could have a CN/HCN ratio as low 
as ~ 3 - 5. In this case the CN/HCN ratio seems to indi- 
cate that IRAS16594-4656 is just starting the post-AGB phase. 
Other PPNe have a CN/HCN ratio up to an order of magnitude 
smaller, e.g., for CRL618, CN/HCN~0.5 and for CRL2688, 
CN/HCN~0.2. IRCh- 10216 has a smaller CN/HCN ratio yet, 
and the more evolved object NGC7027 has a CN/HCN ratio of 
about 10. 

Both IRAS 16594-4656 and IRAS 17 150-3224 have very 
low HCN/CO ratios compared to CRL618, CRL2688 and 
OH231.8-H4.2. IRAS 16594-4656 has an HCN/CO ratio of 
6x10 ^ whereas CRL618 and CRL2688 have ratios of 6x10"^ 
and 3x10-2 respectively (Bujarrabal et al. 1994). IRAS 17 150- 
3224 has an HCN/CO ratio of 2x10"^ which is again smaller 
than the 2x10-^ of OH231.8-H4.2 (Bujarrabal et al. 1994). 

From what can be impUed by Table 8, IRAS 17 150-3224 is 
similar to OH23 1.8-1-4.2 in >^CO, and HCN (despite the spread 
in observed values). The relatively high fractional abundance of 
HCN (half that of the C-rich object IRAS 165 94-465 6) would 
imply that IRAS17150-3224 is at a stage before heavy ioni- 



sation of the molecular matter occurs (and the CN fractional 
abundance rises), i.e., IRAS17150-3224 is at a less evolved 
stage than IRAS 165 94-465 6. 

7.3. Chemical evolution from the AGB 

One has to ask why there seem to be two differing evolution- 
ary paths for PPNe: why are some PPNe molecule-rich, hke 
CRL618, and why are some molecule-poor, like the two objects 
discussed here? Visually, IRAS 16594-4656 and IRAS17150- 
3224 are not very different from CRL618. All three have a sim- 
ilar bipolar structure, with a pinched waist indicating an equa- 
torial density enhancement and possibly a collimating mecha- 
nism. In the case of CRL618, this mechanism could very well 
be a circumstellar torus, the presence of which has been im- 
plied from observations (see Herpin & Cernicharo 2000, and 
references therein). Hence it does not seem too exotic to con- 
sider similar structures for IRAS 165 94-465 6 and IRAS17150- 
3224. Many of the PPNe identified do have such structures, 
which have been directly observed (e.g., IRAS 17 106-3046 - 
Kwok et al. (2000), IRAS 17245-3951 - Hrivnak et al. (1999), 
IRAS04296-I-3429 - Sahai (1999), Hen 401 - Sahai et al. 
(1999), IRAS 1744 1-24 11 - Su et al. (1998); Kwok et al. 
(1996)) and there is strong evidence from polarimetry (Su 
et al. 2000) that a circumstellar disk or torus does exist around 
IRAS 16594-4656. 

CRL618 owes its molecular richness, in a large part, to the 
shielding eff'ects of its torus. This can clearly be seen in the 
chemical model of Woods et al. (2003a) - the high levels of 
radiation are only abated by high initial densities and a slow 
torus expansion velocity, and when the density, or essentially 
the extinction, becomes too low, complex molecules are de- 
stroyed. One obvious implication that this has for PPNe such as 
IRAS 16594-4656 and IRAS 17 150-3224 is that their collimat- 
ing torus must not be dense enough to synthesise and protect 
complex molecules. This brings in an interesting interplay be- 
tween the column density of material required to produce the 
degree of coUimation seen and the column density of material 
required to shield complex molecules from the often very in- 
tense UV fields, cosmic rays and possibly also X-rays. One 
of the prevalent results from Woods et al. (2003a) was that 
incident UV radiation would start to destroy molecules when 
the optical depth reached approximately 10 magnitudes of ex- 
tinction. This equates to a colurmi density of ~1.6 lO^^cm 2. 
For the particular model of CRL618, the density at the inner- 
most part of the torus is ~10^cm"-'. Numerical models (e.g., 
Ignace et al. 1996) show that only an equatorial density en- 
hancement of 2 or more over the polar direction is enough to 
produce asphericity (bipolarity) in a wind-compressed disk sit- 
uation. Typically an average AGB star (e.g., M = 10"^ Mq yr ' 
and Vexp = 15 km s"^) would have a density of around lO' cm"^ 
at a similar radial distance (lO'^ cm). 

To support the hypothesis that these two PPNe are low in 
fractional abundances due to low envelope density, a chem- 
ical model was constructed to follow the evolution of a star 
along the AGB phase and into the PPNe phase. Details of the 
model can be found in Sect. 6.3. The model shows that within 
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~500yr of the cessation of mass-loss UV photons from the 
star start to ionise the molecular material at a very high rate, 
and molecular matter is destroyed rapidly. Evidence of the in- 
creasing UV field can be seen particularly in the profiles of 
CN (Fig. 8), which show an increase in fractional abundance 
at the 500 yr mark, where the UV field starts to become sig- 
nificant. The profiles shown in Fig. 8 are also subtly difi'erent 
from those of steady-state AGB models (compare with the left- 
most set of axes, labelled "Oyr"), indicating that the mildly in- 
creasing UV field does have a marked effect, even in the early 
stages of post-AGB evolution. X-ray photons have little effect 
in this environment, with the chosen assumptions. CO, which 
self-shields in the model, would survive at least 12 000 yr after 
the end of mass-loss at detectable levels. This differs from the 
case of a molecular torus (Woods et al. 2003a) in that in the 
CSE case there is a reasonably slow decline of parent species, 
like HCN, whereas in the torus model, parent molecules are 
destroyed within a hundred years once the extinction provided 
by the circumstellar matter drops to around 10 magnitudes. It 
is also quite marked how rapidly species form in the torus be- 
fore the "radiation catastrophe", in comparison with the slow 
build-up in the CSE model. As a consequence of this, there are 
not the high fractional abundances of photodissociation prod- 
ucts (e.g., CN, C2H) in the CSE case that are found in the torus 
model, and in molecule-rich PPNe, like CRL618. 



Using the fractional abundances calculated from observa- 
tions and the chemical model an approximate post-AGB age for 
IRAS 1 6594-4656 can be estimated. IRAS 1 6594-4656 seems to 
be older than IRAS 1 7 1 50-3224, with quite a high CN fractional 
abundance and also a higher CN/HCN ratio. Through compari- 
son with Fig. 8 in particular, the figure of 370 yr, named by van 
der Veen et al. (1989), seems a good estimate, since a CN/HCN 
ratio of 3-5 (or even 10) is only possible at times less than 
500 yr. IRAS 17 150-3224 is sHghdy more diflicult to qualify, 
since there is no detection of CN. However, this does suggest 
that the CN/HCN ratio is low. The HCN/CO is high (higher 
than IRAS 16594-4656) and hence a post-AGB age of 150- 
210 yr for IRAS 17 150-3224 (Hu et al. 1993; Meixner et al. 
2002) seems likely, with the estimate of van der Veen et al. 
( 1 989), 800 yr, appearing to be much too long to be an accurate 
post-AGB age for this object. 



It is not clear whether IRAS 16594-4656 and IRAS 17 150- 
3224 are older or younger than CRL618 and CRL2688. The 
lack of signs of ionisation suggest that they are younger; the 
reasonably high CN/HCN ratios and low HCN/CO ratios sug- 
gest that they are older. Certainly, IRAS 1 6594-4656 fits into the 
" evolutionary sequence" (of post-AGB objects) of Bachiller 
et al. (1997a) between CRL618 and the young PN, NGC7027. 
The age estimates of van der Veen et al. (1989), Hu et al. (1993) 
and Meixner et al. (2002), above, in comparison with the re- 
sults of Bujarrabal et al. (2001) also tend to favour IRAS 16594- 
4656 and IRAS 17 150-3224 being older, since Bujarrabal et al. 
(2001) give post-AGB ages for CRL618 and CRL2688 as 110 
and 200 yr, respectively. 



8. Conclusions 

Two pro to-planetary nebulae, IRAS 16594-4656 and 
IRAS17150-3224, were observed in a wide range of molecular 
lines, but only detected in a few. Calculating fractional abun- 
dances and upper limits from these fines shows that these two 
sources are molecule-poor in relation to other PPNe, such as 
CRL618, CRL2688 and OH231.8+4.2. As a reason for this 
apparent difference the degree of density of the circumstellar 
torus (or disk) is suggested, with molecule-rich PPNe having 
dense, protective and nurturing tori, and molecule-poor PPNe 
having tenuous, or no, tori. To substantiate this, a model of the 
chemistry in a late-AGB/early-PPN circumstellar envelope, 
with no particular density enhancements (such as a torus), is 
used. As expected, it shows that very few molecules reach 
high fractional abundances in the post-AGB phase, and agrees 
well with fractional abundances calculated from observations. 
This is in stark contrast to the model of Woods et al. (2003a), 
who model a dense circumstellar torus in the PPN phase. 
The usefulness of HCN/CO, HNC/HCN and particularly 
CN/HCN ratios in determining evolutionary age is discussed, 
and using these tools, post-AGB ages for IRAS 165 94-465 6 
and IRAS17150-3224 are given. 
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Appendix A: Molecular line spectra 

Spectra of IRAS 165 94-465 6 and IRAS17150-3224 taken with 
the SEST are presented in Figs. A.l - A.3. 
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Fig. A.l. Molecular line observations of IRAS 16594-4656. 
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Fig. A.2. Molecular line observations of 1RAS16594-4656. 
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Fig. A.3. Molecular line observations of IRAS17150-3224. 



